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Abstract 

Charged lepton-flavor violating processes are extremely rare in the Standard Model, but they are predicted to 
occur in several beyond-the-Standard Model theories, including Supersymmetry or models with leptoquarks or 
compositeness. We present a search for such processes in a sample of 117 x 10 6 T(3S) decays recorded with 
the BaBAR detector. We place upper limits on the branching fractions BF(T(3S) — > e ± r =F ) < 5.0 x 10~ 6 and 
BF(T(3S) — > /i ± r T ) < 4.1 x 10~ 6 at 90% confidence level. These results are used to place lower limits on the mass 
scale of beyond-the-Standard Model physics contributing to lepton-flavor violating decays of the T(3S). 
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I. INTRODUCTION 

In the Standard Model, processes involving charged lepton-flavor violation (CLFV) are unobservable since they are 
suppressed by the ratio {^(ml)/M^ v ) 2 < 1(T 48 0,0. Here A(m£) is the difference between the squared masses of 
neutrinos of different flavor and Myy is the charged weak vector boson mass. Many beyond-the-Standard Model (BSM) 
scenarios predict observable rates for these processes [H, 0| , which may lead to striking experimental signatures and 
provide unambiguous evidence for new physics. There have been considerable efforts both in experimental searches 
for CLFV decays of the \i and r and theoretical predictions for their branching fractions in various BSM scenarios, 
but CLFV in the Y sector remains relatively unexplored. If new particles contributing to CLFV couple to b quarks, 
such processes may be observable in decays of the T. Unitarity-based considerations allow relations to be derived 
between CLFV decay rates of the /i and r and corresponding CLFV decay rates of the Y |5j. In particular, the current 
experimental constraint BF(t — > 11' E') < (2 — 4) x 10~ 8 [1,0], m which I and t! are charged leptons of either the 
same or different flavor, implies BF(Y(3S) — > ^ ± r =F ) < (3 — 6) x 10~ 3 . A decay rate of this magnitude would result 
in O(10 5 ) T(3S) — > ^ ± r =F decays in our dataset and would be easily observable. If the new physics contributing to 
CLFV is in the Higgs sector, it would preferentially couple to heavy quark flavors, further motivating the search for 
CLFV in the bottomonium sector. 

The rates for the CLFV decays T(4S) — > £ ± r T are too small to be observed, since even T(45) — ► t + t~ has not 
yet been observed. However, the branching fractions for rare decays of the narrow Y(3S) resonance are enhanced 
roughly by Tx(4S)/^r(3S) ~ 10 3 with respect to those of the T(45 l ), dramatically increasing the sensitivity to rare 
processes. In this analysis we search for the CLFV decays Y(3S) — > £ ± r =F (£ = e, fj) using data collected with the 
BABAR detector. The prior constraints on CLFV Y decay branching fractions come from the CLEO experiment H, 
which placed the 95% confidence level upper limit BF(Y(3S) -> M ± ^ T ) < 20.3 x ICT 6 . This analysis is more than a 
factor of four more sensitive to this decay and places the first upper limits on BF(Y(3S) — > e ± r =F ). Since the decays 
we are searching for are necessarily mediated by new particles produced off-shell in loops, their measurement probes 
mass scales far exceeding the PEP-II center-of-mass (CM) energy up to the TeV-scale [9j. 

II. THE BABAR DETECTOR AND DATASET 

The data used in this analysis were collected with the BaBAR detector at the PEP-II asymmetric-energy e + e~ 
collider. We search for Y(3S) — > e ± r =F and Y(3S) — > ^ ± r =F decays in a sample of (116.7 ± 1.2) x 10 6 Y(3S) decays 
corresponding to an integrated luminosity of 27.5 fb _1 . Data collected at the Y(4S) corresponding to 77.7 fb _1 and 
data collected 30 MeV below the Y(3S) resonance corresponding to 2.6 fb _1 constitute background control samples 
which are not expected to contain signal events. An additional data sample collected at the Y(3S) corresponding to 
1.2 fb _1 , for which the limit from the CLEO collaboration implies that less than 5 signal events should be present 
per channel, is also used as a background control sample which is not included in the 27.5 fb _1 sample. These control 
samples are analyzed to verify that a signal yield consistent with zero is obtained. 

Simulated events are also produced and analyzed in order to optimize the selection and fitting procedure and to 
compare to data. The background to our events is dominated by continuum QED processes, with an additional 
contribution from resonant Y(3S) — > t + t~ production. The KK2F generator [Tcj is used to produce /i-pair and 
r-pair events while taking into account the effects of initial state radiation. The BHWIDE generator [ll[ is used to 
produce Bhabha events, also taking into account initial state radiation. The EvtGen generator [l^] is used to produce 
generic Y(3S) decays as well as 1.6 x 10 5 signal Y(3S) — > e ± r T and Y(3S) — > /x ± t t decays. The efficiency for qq 
events (q — u, d, s, c) and for two-photon processes to pass selection is found to be negligible, so these processes are not 
included. The simulated /x-pair, r-pair and generic Y(3S) samples correspond to roughly twice the number of events 
in the Y(3S) dataset, while the Bhabha sample, which constitutes a small background to our events, corresponds to 
roughly half the number of events. PHOTOS 13] is used to simulate radiative corrections, and GEANT [14| is used 
to simulate the interactions of particles traversing the BABAR detector. 

The BABAR detector is described in detail elsewhere jl5| . Charged particle tracking is provided by a five-layer double- 
sided silicon vertex tracker (SVT) and a 40-layer drift chamber (DCH). Photons and neutral pions are identified and 
their energy deposition is measured using the electromagnetic calorimeter (EMC), which is comprised of 6580 thallium- 
doped Csl crystals. These systems are mounted inside a superconducting solenoidal coil providing a 1.5-T magnetic 
field. The Instrumented Flux Return (IFR) forms the return yoke of the superconducting coil, instrumented in the 
central barrel region with limited streamer tubes for the identification of muons and the detection of clusters provided 
by neutral hadrons. 



8 



III. ANALYSIS METHOD 
A. SIGNAL SIGNATURES AND BACKGROUNDS 

We search for the decays T(3S) — > e ± r =F and Y(3S) — > ^i ± r =F . The signature for these events are two oppositely- 
charged tracks, a primary electron or muon with CM momentum close to the beam energy, and a secondary charged 
lepton or charged pion from a tau decay. If the tau decays leptonically, we require that the primary and secondary 
leptons are of different flavor. If the tau decays hadronically, we require one or two additional neutral pions from this 
decay. Thus we define four signal channels (here and in the following charge conjugate final states are implied): 



• leptonic er channel: T(3S) - 


-> e ± r =F j r — ► 






• hadronic er channel: T(3S) - 


-> e^r^, t~ — * 


H~ ifi V T j lt~ ifi V T 




• leptonic fir channel: T(3S) 


-> ^ ± t =f , r" -h 


■ erv T v e 




• hadronic fir channel: T(3S) 


->■ ^ ± t t , r- - 






The r(35) -» £M=, r - -> tt" 
channels are omitted due to strong 


-i/ r , r(35) -> 

contamination 


e ± r =F , r~ — > e~v r v e and T(35) — > At ± T =F , r~ — > /u — i/ T 
from Bhabha and /i-pair events. The main source of baclq 


decay 
;round to 



our events comes from r-pair production, which is dominated by continuum production but has a contribution from 
resonant T(3S) — > t + t~ production as well. There is also a background contribution to the T(3S f ) — > e ± r T search 
from Bhabha events in which one of the electrons is misidentified, and to the T(3S) — > /i ± r^ search from /i-pair events 
in which one of the muons is misidentified or decays in flight, or an electron is generated in a material interaction. 

B. STRATEGY AND SELECTION 

Our analysis strategy is to select events with the correct particle types with as high an efficiency as possible and 
then fit for the number of events containing an electron or muon with CM momentum close to the beam energy, which 
are identified as signal events. A preliminary unblinded analysis is performed using the 1.2 fb _1 T(3S) data sample 
in order to ensure agreement between data and simulation and to validate the analysis method. A blinded analysis is 
then performed using the full T(3S) dataset in which events satisfying 0.95 < pi/Eg < 1.00 are excluded, where p\ 
is the electron (muon) CM momentum for the T(3S) — » e ± r T (T(3S) — > /i ± r T ) search and Eg is the beam energy in 
the CM system, equal to y/s/2. The blinding criterion rejects more than 99% of observable Y(3S) -> t±rf decays. 

The event selection proceeds in two steps, a preselection for all four signal channels followed by a channel-specific 
selection. This procedure is designed to eliminate as much of the Bhabha and /i-pair backgrounds as possible while 
retaining high signal efficiency. To pass preselection the event must have two tracks of opposite charge, both consistent 
with originating from the primary interaction point, with opening angle greater than 90° in the CM frame. A neutral 
particle must satisfy Eoep > 50 MeV, where Edep is the energy deposited in the EMC, and must have a transverse 
shower profile consistent with that expected from an electromagnetic shower. The event must satisfy requirements 
designed to select e + e~ — > t+t~ events, which are very efficient for our signal. We require that cos(0^ ss ) <0.9 and 
coa{6^l s ) >-0.9, where 9 l ^ iss (6^^f s ) is the polar angle of the missing momentum vector in the lab (CM) frame. These 
requirements suppress events in which particles are lost because they travel along the beam direction; they also suppress 
Bhabha, /i-pair and two-photon processes. We require that Myis/ ' \fs <0.95, where My is is the mass of the 4- vector 
obtained by adding up the 4-vectors of the two tracks plus those of any additional neutral particles in the event; this 
requirement also suppresses Bhabha and /x-pair backgrounds. Finally, we require (p\ +P2)±/(V~s — \pi \ — \p2\) >0.2, 
where p\ and £>2 are the CM momentum vectors of the two tracks. The requirement on this kinematic variable is 
effective at suppressing two-photon processes as well as suppressing beam-gas interactions, in which a particle from 
the beam undergoes small-angle scattering after interacting with residual gas particles in the beam-pipe. 

Particle identification is performed using a multivariate algorithm [16f | that uses measurements from all of the 
detector sub-systems, including but not limited to E/p and shower profile from the EMC, the specific ionization 
(dE/dx) from the tracking detectors, and the number of hits in the IFR. An electron selector and muon veto, 
combined with the requirement that the particle falls within the angular acceptance of the EMC, are used to identify 
electrons. A muon selector and electron veto are used to identify muons, while a charged pion selector, electron 
veto and muon veto are used to select charged pions. Charged particle misidentification rates are crucial because 
they determine the selection efficiencies for Bhabha and /i-pair backgrounds. The probability for a muon to pass the 
electron selection criteria is O(10~ 5 ). The probability for an electron to pass the muon selection criteria, plus the 
additional requirement that at least 4 IFR hits are associated to the track as required in the leptonic er channel, is 
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TABLE I: Channel-specific selection for the four signal channels. The subscript 1 refers to the primary track which is the 
electron (muon) for the T(3S) — > e ± r T (T(3S) — ■> (i t^) search; the subscript 2 refers to the other track. CM momenta are 
denoted by p, Eb is the beam energy and A(f> CM is the difference between the azimuthal angles of the two tracks in the CM 
frame. The secondary track CM transverse momentum is denoted by px2 , and Nifh2 is the number of IFR hits associated 
to the secondary track. The invariant mass of the vr*^ system is denoted by M(7r ± 7r°). If there are two neutral pions in the 
event, the 7r° giving the ^-k mass closest to m p =0.77 GeV/c 2 is chosen. The requirement on M(-k ± ty°-k°) is included only if 
there are two neutral pions in the event. Empty table entries indicate that no cut is used for the given channel. 



Quantity 


leptonic er 


hadronic er 


leptonic fir 


hadronic /it 


PID 


le, ljj, 


le, 1tt ± , 1 or 2 7r°'s 


le, 1/i 


lfi, 1tt ± , 1 or 2 7T°'S 


Pi/E B 


>0.75 


>0.75 


>0.75 


>0.75 


Pti 1 E B 




>0.05 






P-z/Eb 




<0.8 




<0.8 


A^ CM 






<172° 




NlFR2 


>3 








M(7T ± 7r°) 




0.4-1.1 GeV/c 2 




0.4-1.1 GeV/c 2 


M(7T ± 7r 7r°) 




0.6-1.5 GeV/c 2 




0.6-1.5 GeV/c 2 



TABLE II: Summary of selection for 27.5 fb _1 of T(3S) data and simulated events. Shown are estimated signal efficiencies 
and estimated number of background events, determined from the simulated event samples, and number of data events passing 
selection in the four signal channels. The signal efficiencies include the r branching fractions. All errors are statistical only and 
the errors in the r branching fractions are included in the signal efficiency uncertainties. Systematic effects leading to potential 
discrepancies between the number of data and simulated events passing selection are discussed in Sec. IIVI 



Quantity 


leptonic er 


hadronic er 


leptonic /ir 


hadronic fjbr 


esiG (%) 


4.72±0.05 


4.94±0.06 


4.16±0.05 


6.21±0.06 


Nbkg 


18966±100 


20524±101 


19995±100 


28087±119 


Ndata 


18720 


20548 


19966 


27479 



O(10~ 6 ). The probability for an electron or muon to pass the charged pion selection criteria is O(10 _1 ); this large 
misidentification rate is the reason for requiring additional neutral pions in the event. A pair of photons with invariant 
mass 0.11 GeV/c 2 < M 77 <0.16 GeV/c 2 is selected as a neutral pion candidate. 

The channel-specific selection consists of particle identification requirements, requirements on CM momenta, and 
additional kinematic selection criteria aimed at further suppressing the Bhabha and /i-pair backgrounds as summarized 
in Table [IJ A summary of estimated signal efficiencies, estimated number of background events, and number of data 
events passing selection for each of the four signal channels in the 27.5 fb _1 Y(3S) dataset is displayed in Table HT1 
After including all selection requirements, typical signal efficiencies are 4-6% (including the r branching fraction), and 
typical r-pair background efficiencies are (6 — 10) x 10~ 4 , depending on the signal channel. 

Due to the combination of particle identification selectors and vetoes used in this analysis, it is in general not 
possible for a single event to pass selection for more than one signal channel. The exception is that it is possible for a 
r-pair event decaying to a final state containing an electron and a muon with p e /Es, P^/Eb > 0.75 to pass selection 
for both the leptonic er and leptonic (it signal channels. The probability for an event passing selection for one of the 
leptonic channels to pass selection for the other leptonic channel is determined to be 2%, in which case the event is 
included in both leptonic signal channels. 

C. FIT METHOD 

The discriminant kinematic variable separating signal from background is the beam-energy-normalized primary 
track CM momentum x=p\j Eb, where the primary track is defined to be the electron in the T(3S) — > e ± r =F search 
and the muon in the T(35) -► m ± t =f search. A comparison of the x distributions between selected data and simulated 
events is displayed in Fig.[TJ After selection an unbinned, extended maximum likelihood fit is performed. A probability 
density function (PDF) consisting of a sum of three components is fit to the measured x distribution for each signal 
channel. The three components are signal T(3S) — * i? ± r T production, r-pair background, and Bhabha (//-pair) 
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FIG. 1: The beam-energy-normalized primary track CM momentum distributions x=p\/Es for the selected data and simulated 
events in the four signal channels. Data are denoted by points with statistical error bars, Standard Model background processes 
are denoted by histograms with solid or dotted outlines, and signal events are denoted by histograms without solid outlines, 
where BF{Y(3S) -> ^ ± r T ) has been set to 10" 4 . Solid green is continuum r-pair production and the unshaded region is 
resonant r-pair production. For the er channels, the yellow histogram with a dotted outline is Bhabha events and the solid 
light blue histogram with no outline is signal Y(3S) — > e ± r =F production. For the [it channels, the magenta histogram with a 
dotted outline is /i-pair events and the solid red histogram with no outline is signal T(3S I ) — » fj, production. 



background for the T(3S) — ► e ± r =F (T(3S) — ► /i ± t t ) search, with the yields of the three components left as free 
parameters in the fit. PDFs for signal, Bhabha, and /i-pair backgrounds are determined using x distributions from 
simulated events which are required to pass the same selection criteria as data. The procedure for determining the 
r-pair background PDF is discussed below. The signal yield Nsig is extracted by the fit and is used to measure 
the corresponding branching fraction according to BF = Nsig/( € Sig x ^t(3S))j where esiG is the signal selection 
efficiency and Nx(3S) is the number of produced Y(3S) decays. 

For Y(3S) — > I decays the x distribution is sharply peaked at ir=0.97 because some of the collision energy 
is carried away by the r mass. The x distribution for the signal has a width of about 0.01 and a radiative tail, 
which is more pronounced for the er than for the fiT decays. The resulting x distribution is well-described by a 
modified Crystal Ball function [17], hereafter referred to as a double-sided Crystal Ball function, which has both a 
low-energy and a high-energy tail. The shape of this function is extracted from fits to the simulated Y(3S) — > £ ± r T 
x distributions. The Bhabha and /i-pair backgrounds have a peaking Gaussian component near x=l (approximately 
three times the resolution above the signal peak) also with a width of about 0.01, and a broad component modeled 
by an Argus function [l8| which truncates near x=l. The shape of the Argus plus Gaussian fit function is extracted 
from fits to the x distributions of simulated Bhabha and /i-pair events. 

The PDF describing the r-pair background x distribution is given by the convolution of a polynomial and a detector 
resolution function. The polynomial is given by Eq. [1] 

p{x) = (1 - x/xmAx) + C 2 (l - x/xmAx) 2 +C 3 (1 - xjx M Axf (1) 

in which Xmax determines the kinematic cut-off and C2 and C3 determine the polynomial shape. The detector 
resolution function is extracted from fits to lepton momentum resolution distributions from simulated r-pair events. 
The resulting distributions are found to be well-described by a double-sided Crystal Ball function, whose width varies 
linearly as a function of generated momentum and is extrapolated to 0.96 x Eb, approximately la below the r-pair 
kinematic cut-off. The value of xmax is extracted from a fit to the 77.7 fb^ 1 Y(4S) data control sample as shown 
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LeptonicT(4SHei Channel (x 2 /ndf = 40.8/51) 




Leptonicr(4SHirc Channel (x 2 /ndf = 50.8/51) 




FIG. 2: Fit results for the 77.7 fb" 1 T(4S) dataset for the four 
PDFs, the solid line indicates their sum. The thin green dashed 
line is the Bhahba background PDF for the er channels and the 
a close-up of the region 0.95< x <1.02. 



Hadronic Y(4SHet Channel (x 2 /ndf = 49.5/51) 




Hadronic T(4SHht Channel (% 2 /ndf = 46.4/51) 




signal channels. The dashed lines indicate the two component 
line is the r-pair background PDF, the thick magenta dashed 
p-pair background PDF for the [it channels. The inset shows 



in Fig. [2 while the polynomial shape parameters, which are not strongly correlated with the signal yield, are left as 
free parameters in the fit to T(3S) data. 



D. FIT VALIDATION 



To validate the fit procedure, we perform fits to the T(45) data control sample in order to_verify that a signal yield 
consistent with zero is obtained for each of the four signal channels. The efficiency for B° B° / B + B~ events to pass 
selection is negligible, so that the event compositions for the Y(4S) and T(3S) data samples are similar, with the 
exception that the T(3S) data sample contains an 0(10%) contribution from resonant r-pair production. The Y(4S) 
dataset is divided into a 2/3 (51.8 fb" 1 ) 'control' sample and a 1/3 (25.9 fb _1 ) 'fit' sample, where the size of the fit 
sample is chosen to be comparable to the full T(3S) dataset. The polynomial shape and cut-off are extracted from a 
fit to the 'control' sample, in which the r-pair and Bhabha//x-pair background yields are floated but the signal yield is 
fixed to zero. The fit is then repeated using the 'fit' sample, in which the polynomial shape and cut-off are fixed but 
the yields for signal, r-pair, and Bhabha//i-pair backgrounds are floated. Signal yields consistent with zero within 
±1.4cr are obtained for the four signal channels. We perform similar studies using the 2.6 fb^ 1 T(3S) off-resonance 
sample and the 1.2 fb _1 T(3S) on-resonance sample, for which the limit from the CLEO collaboration implies that 
less than 5 events should be present per channel. These studies confirm the results of the fit validation study using 
the T(AS) data control sample. 

A large number of pseudo-experiments are also performed in order to further validate the fit procedure. Sample 
data sets are generated from the background PDF, simulated signal events are added and the fit is performed using 
the resulting simulated dataset. This procedure is repeated multiple times and the extracted signal yield, error and 



12 



TABLE III: Summary of systematic uncertainties in the signal yield, signal efficiency and number of produced Y(3S) decays. 
The uncertainties in the signal efficiencies and number of Y(3S) decays include both statistical and systematic effects. Also 
shown is the uncertainty due to potential bias in the fit procedure. 



Quantity 


leptonic er 


hadronic er 


leptonic \vr 


hadronic [it 


NsiG 


5.8 events 


8.2 events 


6.7 events 


11.5 events 




2.0% 


3.2% 


2.2% 


2.9% 


Nr(3S) 


1.0% 


1.0% 


1.0% 


1.0% 


Fit Bias 


1.6 events 


1.3 events 


0.8 events 


1.3 events 



pull are recorded. This study confirms that no large biases are observed and that the calculated statistical errors are 
accurate. 



IV. SYSTEMATIC ERROR STUDIES 

The decay branching fractions are determined using the extracted signal yield, estimated signal efficiency, and 
number of produced T(3S) decays. The uncertainties in these quantities are summarized in Table IIIII There is also 
a systematic uncertainty arising from potential bias in the fit procedure. 

The dominant systematic uncertainty in the decay branching fractions comes from the uncertainty in the extracted 
signal yields due to uncertainties in the PDF shape parameters. To assess these uncertainties each parameter pi is 
varied by its uncertainty and the resulting change in signal yield Si—ANsic is determined. The total systematic 

uncertainty in the signal yield is given by Stot = CS, where S =< di • • • 5n > and C is the parameter correlation 
matrix, giving a systematic uncertainty in the signal yield of 6-12 events depending on the signal channel. 

The signal efficiency is determined using simulated events, so there is a systematic uncertainty arising from any 
potential discrepancies between data and simulation. To assess this uncertainty, the relative difference between the 
yields for data and simulated events is taken using r-pair control samples from a portion of the sideband of the x 
distribution defined by 0.8 < x < 0.9. This procedure gives systematic uncertainties in the signal efficiency of 2-3% 
depending on the signal channel, due to uncertainties in the particle identification, tracking, trigger and kinematic 
selection efficiencies. The statistical uncertainty in the signal efficiency is about 1% for each channel. The number 
of produced T(3S) decays in the dataset is estimated to be N r ( 3S ) — (116.7 ± 1.2) x 10 6 by counting the number of 
multihadron events in the dataset, giving a 1% systematic uncertainty. 

To assess the potential bias in the fit procedure, the study described in the previous section using a large number 
of pseudo-experiments is repeated with the generated signal yield fixed to the value extracted by the fit to 27.5 fb _1 
T(3S) data. For those signal channels in which a negative signal yield is extracted, the generated signal yield is 
fixed to zero. The fit is performed multiple times and the deviation from zero in the distribution of extracted minus 
generated signal yield is taken as the uncertainty due to the fit bias, yielding uncertainties of about 1 event count. 

V. RESULTS 

The fit results for the full T(3S) dataset are displayed in Table IIVI and in Fig. [3J After including statistical and 
systematic uncertainties, the extracted signal yields are all consistent with zero within ±2. la. The upper limits, mean 
values, and asymmetric errors on the decay branching fractions are determined by performing a likelihood scan. For 
each of the four signal channels, the signal yield is scanned in steps from -50 to 150, the fit is performed with the 
signal yield fixed at that value, and the likelihood L is extracted. The signal yield is converted to a branching fraction 
according to BF = Nsig/ { e siG x ^r(3S))j gi ym g the likelihood as a function of branching fraction. Systematic 
effects which are uncorrelated between the signal channels are included in the likelihood curves for each channel 
individually, resulting in widening of the likelihood function. The two er (fJ/r) likelihood curves are multiplied to 
obtain the combined er (/ir) likelihood curve, and correlated systematic effects are then included. The branching 
fraction upper limit at 90% confidence level is determined by integrating the likelihood function and finding UL such 
that Jq L L d(BF)/ J °° L d(BF) — 0.9, which is equivalent to a Bayesian upper limit extraction in which the prior is 
taken to be the step function 0(0). The most probable value is taken to be the BF corresponding to the maximum 
of the likelihood function, and the asymmetric errors are determined by finding ±A(BF) such that A(— 2 log(L)) = l. 
The likelihood scan is displayed in Fig.[U and the results of the scan are summarized in Table fVl 
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TABLE IV: Summary of the fit results for the four signal channels of the Y(3S) dataset. Shown are the extracted signal yield 
Nsig, the extracted Bhabha//i-pair background yield Npkbkg, and the extracted r-pair background yield Nbkg- The first 
error is the parabolic statistical error, the second error (displayed only for the signal yield) is the systematic error from the PDF 
shape uncertainties and the fit bias. The signal yields for all channels are consistent with zero within ±2.1ct after including 
statistical and systematic uncertainties. 



Yield 


leptonic er 


hadronic er 


leptonic \xt 


hadronic \it 


Nsig 


21 ±12 ±6 


-1±14±8 


-16±9±7 


42 ± 17 ± 12 


NpKBKG 


25 ± 15 


63 ±16 


35 ±13 


57 ± 12 


Nbkg 


19611 ± 141 


21523 ± 147 


20923 ± 145 


28661 ± 170 



LeptonicT(3SHer Channel fe 2 /ndf = 51 .8/51 ) Hadronic T(3SHet Channel fa 2 /ndf = 40.7/51 ) 




LeptonicT(3SHnt Channel (x 2 /ndf = 35.4/51 ) Hadronic T(3SH>ux Channel (x 2 /ndf = 45.3/51 ) 




FIG. 3: Fit results for the 27.5 fb _1 T(3S) dataset for the four signal channels. The dashed lines indicate the three component 
PDFs, the solid line indicates their sum. The thin green dashed line is the r-pair background PDF, the medium magenta 
dashed line is the Bhahba background PDF for the er channels and the p-pair background PDF for the fir channels, and the 
thick red dashed line is the signal PDF. The inset shows a close-up of the region 0.95< x <1.02. 
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TABLE V: Summary of the results of the likelihood scan. Displayed are the 90% confidence level upper limits and most 
probable values (MPV), with negative and positive asymmetric errors included. 





UL 


MPV 


BF(T(35) -> e ± r T ) (xHT 6 ) 


< 5.0 


9 9+1.9 


BF(T(35) -v fj^r^) (xlO^ 6 ) 


< 4.1 


i 9+1.9 
J -^-1.9 



Combined T(3S) -» ei 




Combined T(3S) -> \vt 




FIG. 4: The results of the likelihood scan. The likelihood curves for the two T(3S I ) — > e ± r T and two T(3S) — > fi signal 
channels have been multiplied to obtain the combined likelihood curves. The red dashed lines indicate statistical uncertainties 
only, the solid blue lines have systematic uncertainties incorporated. The vertical lines bound 90% of the positive integral of 
the total likelihood curve, which is shaded in green. 



VI. CONCLUSIONS 



This paper presents a search for the CLFV decays T(3S) — ► e ± r =F and T(3S) — > /i ± r T . A maximum likelihood 
fit is performed using the primary lepton CM momentum distribution to extract the signal yield. No statistically 
significant signal is observed and the results are used to place the following 90% confidence level upper limits on the 
decay branching fractions: BF(T(3S) -> e ± r = F) < 5.0 x IO" 6 and BF(T(3S) -> m ± t t ) < 4.1 x 10~ 6 . These results 
represent the first upper limits on BF(T(3S) — > e^r^) and represent a sensitivity improvement of more than a factor 
of four with respect to the previous upper limit on BF(T(3S) — ► ^ ± r =F ) Q. 

Effective field theory allows the effects of BSM physics at some large mass scale contributing to CLFV T(3S) decays 
to be parameterized at low energy by a four-fermion interaction with coupling constant ajv and mass scale A. This 
allows the following relation to be derived [lj| [2(| : 



r(T(35) 



r(r(35) -» £+£-) 




,M) (2) 



in which qt, is the charge of the b quark. The dilepton branching fraction of the T(3S) is taken to be the average of the 
PDG [HI values of BF(T(3S) -> e+e~) and BF{T{3S) m + aO, giving BF(Y(3S) -> £+^-) = (2.18±0.15) x IO" 2 . 
Using this result, assuming strong vector coupling with —a^t =1, and taking into account the uncertainty in 
the dilepton branching fractions of the Y(3S), our results place the 90% confidence level lower limits >1.4 TeV 

and A^ >1.5 TeV on the mass scale of BSM physics contributing to CLFV in the bottomonium sector. This result 
improves upon the previous result from the CLEO collaboration [8[ , which used a similar analysis technique to place 
the 95% confidence level lower limit A^ >1.34 TeV. 
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